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Abstract 

Analytic properties of physical quantities in the cosmic fluid such 
as energy density p{t) and Hubble parameter H{t) are investigated 
near the future singularity (Big Rip). Both 4D and 5D cosmologies 
are considered (the Randall-Sundrum II model in the 5D case), and 
the fluid is assumed to possess a bulk viscosity (. We consider both 
Einstein gravity and modified gravity, where in the latter case the 
Lagrangian contains a term i?° with a a constant. If ( is proportional 
to the power (2a — f ) of the scalar expansion, the fluid can pass from 
the quintessence region into the phantom region as a consequence of 
the viscosity. A property worth noticing is that the 4D singularity on 
the brane becomes carried over to the bulk region. 



1 Introduction 



The possibility of crossing the w = —1 barrier in dark energy cosmology has 
recently become a topic of considerable interest. One usually assumes that 
the equation of state for the cosmic fluid can be written in the form 

P = wp, (1) 
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where w is a constant, li w = —1 the fluid is called a "vacuum fluid", 
with peculiar thermodynamic properties such as negative entropy [T]. More 
general forms for the equation of state can be envisaged, such as 

p = w{p)p = -p-f{p), (2) 

which is a form that we shall consider below. As is know, cosmological obser- 
vations indicate that the present universe is accelerating. Recent discussions 
on the actual value of w can be found, for instance, in refs. [2l |3l H]. Perhaps, 
w is even an oscillating function in time. For discussions on time-dependent 
values of w, one may consult Refs. El [7j. The possibility of crossing from 
the quintessence region (—1 < w < —1/3) into the phantom region w < — 1, 
is obviously of physical interest. It may be noted that both quintessence and 
phantom fluids lead to the inequality p + 3p < 0, thus breaking the strong 
energy condition. 

Once being in the phantom region, the cosmic fluid will inevitably be led 
into a future singularity, called the Big Rip [SI El [IHl E] • And this brings 
us to the main theme of the present paper, namely to give an overview of 
the behavior of central physical quantities near the future singularity. This 
is the case of main interest. We think that such an exposition should be 
useful, not least so because the situation is rather complex. Namely, there is 
a variety of different factors at play here: (i) the thermodynamic parameter 
w{p), (ii) the possible time dependence of the bulk viscosity, ( = C(^), and 
(iii) the adoption of Einstein's gravity, or a version of the so-called modified 
gravity. (For an introduction to modified gravity theories, one may consult 
Refs. [Ellis].) 

To begin with, it is convenient to quote from Ref. [TT] the classification 
of possible future singularities: 

(i) Type I ("Big Rip"): For t tg, a — >• cxo, p — *• oo, and \p\ — > oo, or p 
and p are finite aX t = tg- 

(ii) Type II ("sudden"): For t ^ ts, a ^ ag, p — > ps, and \p\ oo, 

(iii) Type III: For t ^ tg, a ^ as, p ^ oo, and \p\ — > oo, 

(iv) Type IV: For t ^ tg, a ^ ag, p ^ 0, \p\ 0, or p and p are finite. 
Higher order derivatives of H diverge. 



Here the notation is standard, a meaning the scale factor and tg referring 
to the instant of the singularity. The above classification was introduced in 
the context of ideal, i.e., nonviscous, cosmology. We can however make use 
of the same classification also in the viscous case. 
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In the following, we will present salient features of 4D, respective 5D, 
viscous cosmology theory, and thereafter focus on the classification of the 
various alternatives. 



2 Viscous 4D theory: Basics 

We include this basic material mainly for reference purposes. We consider 
the standard FRW metric, 

ds"^ = -dt'^ + a^{t)dx^, (3) 

and set the spatial curvature k, as well as the 4D cosmological constant A4, 
equal to zero. The Hubble parameter is H = d/a, the scalar expansion is 
= U^.fj, = 3H with the four- velocity of the fluid, and k| = 87rG4 is the 
gravitational coupling. Of main interest are the {tt) and (rr) components of 
the Friedmann equations. They are 

e' = 34p, (4) 

'^ + H^ = -^lp^ (5) 

where p = p — (6 is the effective pressure. From the differential equation for 
energy, T^'^.^, = 0, we get 

P+{p + p)e = C9'. (6) 

As a consequence of positive entropy change in an irreversible process, we 
must require the value of C to be non-negative. From the equations above 
we obtain the following differential equation for the scalar expansion, 

0-1^0^-^^,1(^0 = 0. (7) 

In view of the relationship 6 = (y/S K4/2)p/ ^/p we can alternatively reformu- 
late this equation as an equation for the density. 



p-K,^3pf{p)-3KlCp = 0. (8) 
The solution is (cf. Eq. (9) in [T4] ) 

(9) 



73^:4 JP* V^/(P)[1 + '^4CV^//(P)]' 
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We here let t = be the initial (present) time, and let the correspond- 
ing initial density be p*. The functional form of the bulk viscosity ( is so 
far unspecified. The shear viscosity is omitted, due to the assumed spatial 
isotropy in the cosmic fluid. Note the dimensions: [^4] = cm^, [/(p)] = 
[p] = cm~^, [(] = cm~^. Viscous cosmology are treated at various places, for 
instance, in Refs. [SlIiailSlinillSlIISlEniEIlElESlISaES]. 



3 Specific cases in 4D 

We are now in a position to discuss various cases in 4D explicitly. We have 
to distinguish between several alternatives: (i) use of Einstein or modified 
gravity; (ii) possible density dependence of the thermodynamic parameter 
w = w{p); and (iii) possible time dependence of the bulk viscosity ({t). 

3.1 Einstein gravity, w and ( being constants 

Let /(p) = ap, with a a constant. The equation of state is then 

p = wp = —(1 + a)p. (10) 
We can now solve explicitly for the Hubble parameter [131 [26] , 



17 pt/tc 

where H^, is the present-time value of H and tc is the 'viscosity time', 

^-^C)". (12) 

From Eq. (fTTl) it is seen that H{t) becomes singular when the denominator 
vanishes. Let us first for reference purposes set C = 0: 

The nonviscous case. If tgo designates the singularity time, we have 



Then [26], 



(13) 



H{t) = ^ (14) 
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a{t) 



Pit) 



sO 



(15) 



(16) 



The viscous case. If now t^c denotes the singularity time, we get from 
Eq. (dH) 

2 1 

tJn 1 



3q; 



and 



HA 



tsO ~ t 



t 



(17) 



(18) 



Close to the singularity we thus obtain the same singular behavior as in the 
nonviscous case. Moreover, we get the following forms, 

a(t) ~ (t,c-t) t^tsc, (19) 
p{t) ~ (t,^ - t)-\ t tsc- (20) 
The viscosity tends to shorten the singularity time, 

ts( < tso, (21) 

but it does not modify the exponents in the singularity. The singularity is of 
Type I if a > 0, and of Type II if a < 0. 



3.2 Einstein gravity, f{p) = Ap^, and ( being constant 

We shall assume that (3 > 1. From Eq. (Q it is apparent that the last term 
in the denominator dominates for large p. Near the singularity we obtain the 
form 

p(t) ~ (t,C-i)^, t^tso (22) 

which generalizes Eq. (1201) and reduces to it when (3=1. Thus p ^ oo 
implying, according to Eq. ([2]), that also \p\ oo. The Hubble parameter 
becomes 

Hit)r^it,^-t)w^K (23) 

If (3 > 1, a as [a finite value) when t —>■ t^c- The singularity is of Type 
III. If /5 = 1, the singularity is of Type I. 

The material of this subsection was discussed also in Ref. [26], whereas 
the equation of state corresponding to fl22|) and fl23|) was discussed by Nojiri 
and Odintsov [27]. 
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3.3 Modified gravity, w being constant, and ( = t9 

Consider now the following gravity model, 

S = ^l (fx^UoR"" + L^). (24) 

where /o and a are constants, being the matter Lagrangian. This model 
has been considered before; cf., for instance, Refs. [2H1 [23 EOl E3] • The case 
/o = 1 and a = 1 yields Einstein's gravity. The equations of motion following 
from the action above are 

-^/o^?^.^" + «/oi?/..i?"^'-«/oV^V./?"-^ + «/o^^.V'i?"-^ = kIT^,, (25) 

where T^y corresponds to the term in the Lagrangian. For the cosmic 
fluid we have 

Tf,u = pU^JJy + p/i^^, (26) 

where h^^ = g^^ + U^Uy is the projection tensor and p = p — (9 the effective 
pressure. In comoving coordinates, = 1, = 0. We assume now the 
simple equation of state given in Eq. ([1]). 

Of main interest is the (OO)-component of Eq. (l25l) . Using that R = 
6{H + 2H^), Too = P, as well as the energy conservation equation ([6]) which 
in turn follows from V^T^jy = 0, we obtain 

^7/oi?" + 3afo[2H - 3j{H + H^)]R''~^ + 3a(a - l)/o[(37 - l)HR + R]R''-^ 

+ 3a(a - l)(a - 2)foR^R'^-^ = 9kICH, (27) 

with 7 = w + 1. The important point now is that this complicated equation 
for H{t) is satisfied with the following form 

H = H^/X, where X = 1- BHJ, (28) 

B being a nondimensional parameter. For Big Rip to occur, B has to be 
positive. 

Taking the bulk viscosity to have the form 

( = t9^^-^ = r(3if)2"-^ (29) 

with r a positive constant, the time- dependent factors in Eq. fl27j) drop out. 
There remains an algebraic equation, determining B. 
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Of main interest is the time- dependent forms 

C = r{?>HjXf^-\ p = p,/X2° (30) 

As an example, the case a = 2 turns out to yield a cubic equation for B. 
There is one positive root (assuming /o positive), leading to a viscosity- 
generated Big Rip. If a < 0, typically a = —1, there may still be positive 
solutions for B implying that H = H^/X is diverging. By contrast, ( oc 
^-{2a~i) ^j^^ p ^ go to zero. 



4 Relationship to 5D viscous theory 



Let us investigate the possible link between the 4D theory above and the 
analogous viscous theory in 5D space. To this end we consider a spatially 
flat {k = 0) brane located at the fifth dimension y = 0, surrounded by an 
anti-de-Sitter (AdS) space. If the 5D cosmological constant, called A, is 
negative, the configuration is that of the Randall-Sundrum II model (RSII) 
[31] . The 5D coordinates are denoted = (t, x, y), and the 5D gravitational 
coupling is = SnG^. The Einstein equations are 



Rab — -QabR + Qab^ = f^l Tab, 



(31) 



and the metric is 



(32) 



where n{t, y) and a{t, y) are to be determined from the Einstein equations. 

Of main interest are the (tt) and (yy) components of the field equations. 
They are 




1^1 Ttt, 



'^s -^yy 



(33) 



(34) 



(cf. for instance, Refs. [321 [331 [311 [191 [26] ) . Overdots and primes mean deriva- 
tives with respect to t and y respectively. On the brane y = we assume 
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there is a constant tension a, and an isotropic fluid with time- dependent 
energy density p = p{t). The energy-momentum tensor is now 



Tab = 5{y){-a6f,^ + pU^U^ + phf,^)5'^5%. 



(35) 



Applying the junction conditions across the brane we obtain, for arbitrary y, 
after integration with respect to y [35j , 



-)■ 

naj 



1 



C 



(36) 



(37) 



We here let subscript zero refer to the brane. On the brane, no(t) = 1. Recall 
that A and a are constants, and that Eq. fl37j) is a 5D, not a 4D, equation. 
Its essential new feature is that it contains a term. The equation functions 
bridge between 4D and 5D cosmologies. 
We observe the solution for ao(t) if p = 0: 



ao(t) 



X 



6 36 ^ ' 



(38) 



normalized such that ao(0) = 1. 

Inserting p = p*/X^" into Eq. fl37|) we get 
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-A 



'3. 

36 



p* 



1 - BHJ) 



2a 



(39) 



Near the Big Rip, tg = 1/ [BH^), the quantities A and a become unimportant 
and we get 

(/t^/6)p. 



ao(t) ~ exp 



(40) 



(2a- l)(M,)2-(t,-t)2"-i 

showing that if a > 1/2, aQ{t) has an essential singularity. Einstein's gravity 
corresponds to a = 1. The singularity becomes stronger, the higher is the 
value of a. If a < 1/2, ao(t) does not diverge at tg. From Eq. (136!) . 



6A 



6A 



cosh(2py) 



3p 



sinh(2p|?/| 



(41) 



8 



with fi = y — A/6. The important point here is that the Big Rip divergence 
on the brane becomes transferred to the bulk. The bulk scale factor a{t, y) 
diverges for arbitrary y ai t = tg ii ao{t) diverges at ts. There is no funda- 
mental difference between an Einstein fluid and a modified gravity fluid in 
this respect; their behavior is essentially the same. 

In summary, we have discussed viscous dark energy as a particular rep- 
resentative of inhomogeneous equation-of-state fluids and the appearance of 
finite-time future singularities for such energies. It is of interest to note 
that due to the relationship between modified gravity and inhomogeneous 
equation-of-state ideal fluids [36], our findings may be useful in the study of 
future singularities in modified gravity [37] . 
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